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Abstract. We describe a method for bounding the rank of an ehiptic 
curve under the assumptions of the Birch and Swinnerton-Dyer conjecture 
and the generahzed Riemann hypothesis. As an example, we compute, 
under these conjectures, exact upper bounds for curves which are known to 
have rank at least as large as 20, 21, 22, 23, and 24. For the known curve of 
rank at least 28, we get a bound of 30. 



1. Introduction 

Determining the rank of an elliptic curve is a difficult problem, and there 
is currently no known unconditional algorithm for determining the rank of a 
given curve. The basic method for rigorously determining the rank of a curve 
is to find an upper bound for the rank by computing the size of some Selmer 
groups and to find a lower bound for the rank by finding enough independent 
rational points. In theory, if one continues this process long enough, and the 
Shafarevich-Tate group of the curve is finite, the upper and lower bounds 
should eventually coincide and the rank will be determined exactly. 

In practice, things are not so simple. Finding points on the curve is sometimes 
not too bad, but the upper bounds for the rank are more problematic. Even 
the computation of the 2-Selmer rank is difficult, and it becomes prohibitively 
time consuming as the coefficients of the elliptic curve grow; it is easy to write 
down a curve for which the state of the art program for computing the 2-Selmer 
group, John Cremona's mwrank |4j, will effectively take "forever." 

If one is willing to accept the Birch and Swinnerton-Dyer conjecture that the 
rank of an elliptic curve is the same as the order of vanishing of its L-function 
at the central point, then it is possible to use the L-function to get information 
about the rank. In fact, when the order of vanishing is between and 3, it can 
be possible to compute the L-function to enough precision and use some extra 
information about the curve to determine the analytic rank exactly, as is done 
in [2], for example. When the rank is larger than this, though, currently the 
best one can do is determine that the first r derivatives of the L-function are 
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very close to 0, and the (r + l)-st is not, which will provide a very good guess 
for the rank and a rigorous upper bound, assuming BSD. 

This approach has its own problems, as it is much easier to write down a curve 
of large conductor than it is to compute the L-function of such a curve. For 
example, the known curve of rank at least 28 [6], which we will write down 
later, has conductor ^ 3.5 x 10^"^^, and current methods (such as those 
described in [16]) typically require summing on the order of \/N terms to 
compute the central value of the L-function. (It would take a compute about 
10^^ cpu-years just to add 1 to itself 10™ times.) 

We present here a third method which is rather effective at bounding the 
rank, especially when the rank is large compared to the conductor, as long 
as one is willing to assume both the Birch and Swinnerton-Dyer conjecture 
and the Riemann Hypothesis for the L-function of the curve. This method is 
not completely new. It is based on Mestre's method [13] for (conditionally) 
bounding the rank of an elliptic curve based only on its conductor, and it 
was used by Fermigier [8] to study ranks of elliptic curves in certain families. 
However, it does not seem to have gained much traction and does not seem to 
have been used much, if at all, since. 

The idea, in brief, is as follows. Take f{x) to be a function such that /(O) = 1 
and f{x) > for all real x. Then, assuming the Riemann hypothesis, the 
sum fil)-! where 1/2 + ^7 runs over the nontrivial zeros of L(s, L'), will be 
an upper bound for the analytic rank of E. Moreover, for certain choices of 
f{x) this sum may be efficiently evaluated using the explicit formula for the 
L-function attached to E. 

This method has recently been implemented by the author, and is available 
as part of William Stein's PSAGE [18] add-ons to Sage [19]. As an example, 
of what it can do, we will examine 6 curves that are known to have rather 
large rank. We denote these curves as E^, n = 20, 21, 22, 23, 24, 28, where n 
is a known lower bound for the rank. We will write down these curves later 
(they are all taken from A. Dujella's website [5], and at the time of discovery 
each held the record for the curve with largest number of known independent 
rational points). The exact rank is not known for any of these curves. However, 
conditionally we may claim 

Theorem 1.1. Assuming BSD and GRH, En has rank exactly n for n = 
20,21,22,23, and 24, while E28 has rank 28 or 30. 

1.1. Acknowledgements. Most of the computations in this paper run in a 
short amount of time, and were done on the author's personal computer. Some 
longer computations were run on the sage cluster at the University of Wash- 
ington, supported by NSF grant DMS-0821725, and the riemann cluster at the 
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University of Waterloo, funded by the Canada Foundation for Innovation, the 
Ontario Innovation Trust, and SGI. 

The source code for our implementation is available as part of PSAGE [18]. It 
uses Sage [IH], and hence PARI [22], to compute for bad primes, and uses 
Andrew Sutherland's smalljac [20] to compute all other values of a^. 

Parts of this work began while the author was in residence at the Mathematical 
Sciences Research Institute during the Arithmetic Statistics program. Spring 
2011, during which time the author was partially supported by NSF grant 
DMS-0441170, administered by MSRI. Discussions during the informal "ex- 
plicit formula seminar," especially with David Farmer and Michael Rubinstein, 
were influential in encouraging this work. 

Currently the author is supported by NSF grant DMS-0757627, administered 
by the American Institute of Mathematics. 

2. Bounding ranks 

2.1. The method. Let 

oo 

n=l p 

be the L-function of an elliptic curve, normalized so that the completed L- 
function A(s, E) = eA(l — s,E), and let c„ be defined by 

L(s,E) ~ ^ n'' 

^ ' 71=1 

More explicitly, if we define and by 

= (1 - a(pK^)(l - 

(note that a and /3 are only well defined up to permutation, and that at least 
one of them will be when p is a prime of bad reduction), then 

(1) c,™ = («(pr+/3(pr)iogp, 

and c„ = when n is not a prime power. 

Our main tool will be the explicit formula for L{s,E), which we state in a 
friendly form in the following lemma. 

Lemma 2.1. Suppose that f(z) is an entire function with f{x + iy) <ti x^^^'^^'^ 
for \y\ < 1 + e, for some e > 0, and that the Fourier transform of f 

POD 

fiy) = / /(x)e-2™^dx 
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exists and is such that 

'logn 



E Cn p /logr; 

n=l ^ 



converges absolutely. Then 



n=l 



where 1/2 + ^7 rans over i/ie nontrivial zeros of L{s,E), where E is an elliptic 
curve with conductor N . 

Proof. A proof of the explicit formula in this form, or in a similar form, can 
be found in various sources, e.g. [101 Theorem 5.12], so we give only a brief 
sketch. The idea is to integrate the function 

L{s,E) 

where -F(l/2 + is) = f{s), on a vertical line to the right of the critical strip 
and, in the reverse direction, on a vertical line to the left of the critical strip. 
By the residue theorem, this integral will be equal to 271^^/(7). One now 
applies the functional equation to write the integral in the left half-plane as 
an integral in the right half-plane. 

The sum over the Fourier coefficients of / arises from shifting contours to the 
region of absolute convergence and using the Dirichlet series for L'{s)/L{s), 
while the other terms arise from shifting the remaining integrals to the line 
3fJ(s) = 1/2. 

The conditions on f{z) are exactly those needed to make sure that this process 
can go through without trouble. Of course, it is also important that L{s, E) 
is entire and that it satisfies a functional equation |23l [211 [I]- ^ 

A convenient function to use in an application of the explicit formula is 

'sin(A7rz) 



f{z) = fiz;A) = 
which has the simple Fourier transform 



Anz 



/(x;A)=(i) (1- 



Ixl < A. 
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With this choice of /, equation ([2]) takes the form 

p<exp(27rA) fc=l ^ ^ ^ 

Since /(7; A) > as long as 7 is real, and /(O; A) = 1, equation ([3]) will give 
an upper bound for the order of vanishing of L{s,E) at s = 1/2, as long as 
the Riemann Hypothesis holds for L{s, E). And if A is not too large, we can 
quickly evaluate the right hand side of equation ([3]) to calculate this upper 
bound. It is also worth noting that, assuming RH, 

-il-^ E E irMpf + B(pf)[i--l 

p<exp(27rA) k=l ^ ^ 

= Olds=l/2L{s,E) 

SO that, in principle, we should be able to get as good a bound for the rank as 
we like through this method. However, as the length of the prime sum grows 
exponentially in A, this method quickly becomes infeasible once A gets a little 
large than 4. 

2.2. Some curves. As an example, we examine 6 elliptic curves from Dujella's 
online tables. They are 

£^20 ■.y^ + xy = x^- 431092980766333677958362095891166X 

+ 5156283555366643659035652799871176909391533088196, 

E2i:y^ + xy + y = x^ + x^ - 215843772422443922015169952702159835x 
- 19474361277787151947255961435459054151501792241320535, 

E22:y^ + xy + y = x^ - 940299517776391362903023121165864a; 

+ 10707363070719743033425295515449274534651125011362, 

E23:y^ + xy + y = x^- 19252966408674012828065964616418441723x 

+ 32685500727716376257923347071452044295907443056345614006, 



E2i:y^ + xy + y = x^- 120039822036992245303534619191166796374x 
+ 504224992484910670010801799168082726759443756222911415116, 
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Curve 


hgNE 


A 


E./(7;A) 


2ttA 


E20 


170.09 


2.0 


21.70 


13.54 


E21 


196.68 


2.5 


22.68 


12.52 


E22 


182.72 


2.0 


23.71 


14.54 


E23 


205.06 


2.5 


24.49 


13.05 


E24 


219.93 


2.5 


25.57 


14.00 


E28 


325.90 


3.2 


31.30 


16.21 



Table 1 . Computed upper bounds for the ranks of some curves, 
along with a heuristic guess of what these bounds should for a 
typical elliptic curve. The sum over the zeros here is rounded 
up; other numbers are rounded to nearest. 



and 



^2 s 2 /20067762415575526585033208 x 10^^\ 

E28 '. V+xy + v = x— X— [ ]x 
^ ^ ^ V+ 2093385427509302303121789565027 

3448161179503055646703298569039072037485594 x 10^° 

+ 4359319180361266008296291939448732243429 



+ 



Each En has n known independent rational points of infinite order, so thus has 
at least rank n. (See [El [111 [U [TTl [121 [6] , or Wi quick reference.) Using the 
methods described above, we compute rank bounds for each of these curves. 
These are listed in Table [H The global root number can be computed for each 
curve. (In Sage, E.root_number(), which uses PARI [22j, will finish quickly 
for £'20, E21, and £'22 and within a few hours for £^23 and E24. For E28 it 
is best to see the mailing list discussion which gives the factorization of the 
discriminant j^.) In each case the root number agrees with the parity of the 
known number of independent points, so to get a tight upper bound for the 
rank we only need to get within 2 of the number of known independent points, 
and so the computation in Table [1] gives the proof of Theorem 11.11 



2.3. Curves of small conductor. For further testing, this method was also 
run on all elliptic curves up with conductor below 180000 (from Cremona's 
tables [3]) using A = 2.0, a computation which ran in under a day on a fast 8 
core computer. In this range there are 790677 isogeny classes of elliptic curves. 



and for all but 9882 isogeny classes it turns out that 



V /(7;2.0) =mnk{Ey, 



in the remaining cases. 



^ /(7; 2.0) = rank(£') + 1, so consideration of the 



root number of the curve gives the exact rank. 
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3.1. Some evidence towards BSD. There is a way in which these com- 
putations can be seen as giving mild evidence in support of the Birch and 
Swinnerton-Dyer conjecture. The upper bound computed for a curve E is the 
value of the sum /(7; A), and as /(7; A) decays fairly rapidly as 7 grows, 
one does not expect this sum to be very large for a typical elliptic curve. 

To obtain a crude approximation to what we might expect the value of this 
sum to be, consider that the local zero density of a typical L{s, E) near the 
central point is approximately . Then, if the zeros are spaced uniformly 
at random (an assumption that is not really correct, but is close enough to 
true for our crude purposes), we might expect that 



possibly with a small adjustment to take into account the parity of the rank. 
(More precisely, we might expect that if we average this sum over all elliptic 
curves of conductor close to iV^, the answer will not be too far from this 
integral.) Thus, when this sum is significantly larger than this estimate, it 
indicates an extreme concentration of zeros near the central point. (It is also 
possible to arrive at more refined version of this heuristic by considering the 
explicit formula. In such it is necessary to assume that the family of 

elliptic curves considered is large enough that ap{E) averages to zero for each 
p, and we notice that the integral of the F-factor plays a small role as well.) 

As some further small evidence for this heuristic, we note that the average of 



over all isogeny classes up to 180000 is approximately .9638. The small dif- 
ference from 1 should be accounted for by the F-factor, which tends to push 
zeros away from the central point. 

It should also be possible to refine this heuristic somewhat to make a guess as 
to what the sum should be for a high rank curve by making the assumption 
that a zero of high order at the central point will push other zeros away. 

3.2. Correctness tests. The method described here is simple enough that it 
is easy to implement, which reduces the likeliness of bugs. It is still important 
to test it where possible, however, in order to have more confidence in its 
correctness. 



3. Further comments 
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As described in Section I2.3[ this code was run on every isogeny class up to 
conductor 180000, and the results there suggest a high degree of confidence in 
the results elsewhere. As a further test, one can also compute many zeros for 
the L-function of an elliptic curve of small conductor, compute the sum over 
zeros directly, and verify that it agrees with our explicit formula implementa- 
tion. This was done with the elliptic curve "llal" for a few values of A, and 
little over 200000 zeros (computed using M. Rubinstein's Icalc package |T7]). 
and the agreement is generally to within about 10~^, which is in fine with 
what is expected using only 200000 zeros, and which is roughly the precision 
to which the integral in the explicit formula was calculated. Similar tests have 
also been done with a smaller number of zeros for other L-functions. 
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